version February 27, 2013: fm 

I Searching for 7-ray blazar candidates among the unidentified INTEGRAL 
"(S|"' 2 sources 

^ ; 3 F. Massaro\ A. Paggi^, R. D'Abrusco^ & G. Tosti^-'^. 

$H ! 

c3 . 4 SLAC National Laboratory and Kavli Institute for Particle Astrophysics and Cosmology, 2575 
^ ! 5 Sand Hill Road, Menlo Park, CA 94025 

o\ '. 

,-H . 6 Harvard - Smithsonian Astrophysical Observatory, 60 Garden Street, Cambridge, MA 02138 
[j^ '- 7 Dipartimento di Fisica, Universitd degli Studi di Perugia, 06123 Perugia, Italy 

^ • 8 Istituto Nazionale di Fisica Nucleare, Sezione di Perugia, 06123 Perugia, Italy 

^ ■ 

O 9 ABSTRACT 

^ ■ 

< 10 The identification of low-energy counterparts for 7-ray sources is one of the 

^ ■ biggest challenge in modem 7-ray astronomy. Recently, we developed and success- 

>■ ■ fully applied a new association method to recognize 7-ray blazar candidates that could 

^^ • be possible counterparts for the unidentified 7-ray sources above 100 MeV in the sec- 

^ . ond Fermi Large Area Telescope (LAT) catalog (2FGL). This method is based on the 

• . Infrared (IR) colors of the recent Wide-Field Infrared Survey Explorer (WISE) all-sky 

O ■ survey. In this letter we applied our new association method to the case of unidentified 

^ . INTEGRAI sources (UISs) listed in the fourth soft gamma-ray source catalog (4IC). 

L; . Only 86 UISs out of the 113 can be analyzed, due to the sky coverage of the WISE 

'Jp . Preliminary data release. Among these 86 UISs, we found that 18 appear to have a 

^ . 7-ray blazar candidate within their positional error region. Finally, we analyzed the 

Swift archival data available for 10 out these 18 7-ray blazar candidates, and we found 

that 7 out of 10 are clearly detected in soft X-rays and/or in the optical-ultraviolet 

band. We cannot confirm the associations between the UISs and the selected 7-ray 

blazar candidates due to the discrepancies between the INTEGRAI and the soft X-ray 

spectra. However, the discovery of the soft X-ray counterparts for the selected 7-ray 

blazar candidates adds an important clue to help understand their origin and to confirm 

their blazar nature. 

II Subject headings: galaxies: active - galaxies: BL Lacertae objects - radiation mecha- 
12 nisms: non-thermal 
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One of the main scientific objectives of the INTEGRAL mission is performing a survey of the 
sky in a mostly unexplored region of the electromagnetic spectrum: from the hard X-ray to the soft 
7-ray band. Since its launch in October 2002, INTEGRAL ha s used the unprecede nted imaging 
capabilities of IBIS (Imager on Board INTEGRAL Sp acecraft: lUbertini et al. 2003b coupled with 
those of ISGRI (INTEGRAL Soft Gamma-Ray Imager: Eebrun et al. 20031) . Combining data from 
these two instruments, it is possible to generate images of the sky with a 12 arcmin (Full Width 
Half Maximum, FWHM) resolution with typical source location accuracy of ~ 1-3 arcmin over a 
~19 degree (FWHM) field of view in the energy range 17-1000 keV. 

The fourth soft 7-ray source catalogll (|Bird et al. 2010|) (hereinafter 4IC) obtained with the 
IBIS 7-ray imager on board the INTEGRAL satellite lists 723 hard X-ray/soft 7-ray sources. In par- 
ticular, the 4IC substantially increased the extragalactic sky coverage including b oth transients an d 
faint persistent objects that can only be revealed with long exposure observations (JBirdet al. 2010|) . 



26 Several observations at low energies have been alr eady performed to decrease the number 

27 of the unidentified INTEGRAL objects (UISs) (see e.g., Masetti et al. 20081 : Masetti et al. 2009 : 



Masetti et al. 2010l . for optical and X-ray observations of UISs); however, a considerable fraction 



of the 4IC sources are still completely unidentified. According to the 4IC, there are 113 UISs, cor- 
responding to about 16% of the whole catalog, and 178 other sources have uncertain classification. 

31 The largest fraction (i.e., 35%) of the associated INTEGR AL sources are A ctive Galactic Nuclei 

32 (AGN), compared to 31% identified as Galactic sources ([Birdet al. 2010|) . For comparison, the 

33 58-month catalog of observations with the BAT hard X-ray detectorEl on board the Swift observa- 

34 tory, lists 109 2 objects detected in the 14-195 keV energy rang e, with 86 unidentified hard X-ray 

35 sources listed (jCusumano et al. 2010l:lBaumgartner et al. 20101) . 



36 Recently, using the WISE all-sky IR survey, we discovered that blazars, the la rgest known 7- 

37 ray class of AGN, can be separated from other extragalactic sources using IR colors ([Massaro et al. 2011a , 

38 hereinafter Paper I). We used the magnitudes of the recent WISE IR all-sky survey performed at 
3.4, 4.6, 12, and 22 /im with an angular resolution of 6.1, 6.4, 6.5 & 12.0 arcsec and with 5a 
point source sensitivities achieving 0.08, 0.1 1, 1 and 6 mJy, in unconfused regions on the ecliptic, 
respectively. The absolute (radial) differences between WISE source-peaks and "true" astrometric 
positions anywhere on th e sky are no large r than ~ 0.50, 0.26, 0.26, and 1.4 arcsec in the four 
WISE bands, respectively dCutri et al. 201lli pl. 



'http://irfu.cea.fr/Sap/IGR-Sources/ 
^http://heasarc. nasa.gov/docs/swift/results/bs58mon/ 



^http://wise2.ipac. caltech.edu/docs/release/prelim/expsup/sec2Jg.html 



44 Moreover, we investigated a sample of blazars d etected by WISE and a ssociated with Fermi- 

45 LAT sources to characterize their IR-7-ray properties ([D'Abrusco et al. 20121 hereinafter Paper II). 

46 This was the first step to develop a new association method for the unidentified 7- ray sources, able 

47 to recognize if there is a 7-ray blazar candidate within their positional error region (|Massaro et al. 2012a 



Massaro et al. 2012bl hereinafter Paper III and Paper IV, respectively). With this new IR diagnos- 



49 tic tool, we searched for 7-ray blazar candidates within the unidentified 7-ray source sample of the 

50 2FGL 7-ray catalog, and for the first time we have been able to provide a candidate counterpart 

51 for 187 out of 313 unidentified 7-ray sources analyzed, having the same IR properties as the 7-ray 

52 emitting blazars (see Paper IV). 

53 In this letter, we apply this association procedure to test whether there is a possible 7-ray 

54 blazar candidate for the UISs using their WISE IR colors. For the selected 7-ray blazar candidates, 

55 we also search Swift pointed observations for the presence of an optical-UV and/ or X-ray counter- 

56 part. This letter is organized as follows: in Section |2] we describe the UIS sample selected for our 

57 investigation; in Section [3] we illustrate the basic details of our new association method, that, in 

58 Section HI we apply to the UISs. Section [5] is devoted to the optical-UV-X-ray counterparts in the 

59 Swift observations available. Our results are discussed in Section [6l 



60 



2. The sample selection 



61 In th e 4IC there are 1 13 sources that are completely unidentified (i.e.. Type =?, according to 

62 Table 3 in lBirdet al. 20101) . while there a re 97 sources tha t are indicated as unidentified transients 

63 (i.e.. Type =?,T, according to Table 3 in iBirdetal. 20101). In addition, w ithin the 4IC, there are 

64 also 32 AGN sources of uncertain type (i.e.. Type =AGN? iBirdet al. 2010l) and another 49 objects 

65 with uncertain classification. 

66 In this letter, we only considered the 86 UISs out of 1 13 listed in the 4IC that lie in the portion 

67 of the sky covered by the WISE Preliminary Source catalog. 
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3. The WGS association method 



In Paper III, us i ng a subsample of the ROMA-BZCAT blazar catalog ([Massaro et al. 2009 : 



Massaro et al. 2010t iMassaro et al. 201 lb|) . detected by WISE and associated with Fermz -LAT 



71 sources (Papers II and III), we presented the parametrization of the WISE 7-ray strip {WGS) based 

72 on the strip parameter s. This parameter, ranging between 0-1, provides an estimate of the dis- 

73 tance between the WGS and the location of a generic WISE source in the IR color parameter space, 

74 and it is weighted for the errors on all the IR colors. We distinguished between WISE sources that 



75 lie in the subregion of the WGS occupied by the BZBs and BZQs using the Sb and Sg parameters 

76 separately (Paper III). 

77 In Paper IV, we presented the association method based on the WGS parametrization. For each 

78 unidentified 7-ray source we defined the searching region corresponding to a circular region of 

79 radius R equal to the semi-major axis of the elliptical source loc ation region at 99.9 99% confidence 



level, centered on the 7-ray position given in the 2FGL catalog (|Nolan et al. 2012h 



81 We calculate the IR colors for every WISE source that lies within the searching region as well 

82 as their Sb and Sg parameters. Given the distributions of generic WISE sources in random regions 

83 of the sky, we distinguish three classes of 7-ray blazar candidates on the basis of their Sb and/or Sg 

84 values: 

85 • class A: WISE sources with 0.24 < Sb < 1.00 and 0.38 < Sq < 1.00; 

86 • class B: WISE sources with 0.24 < Sb < 1.00 or 0.38 < Sg < 1.00; 

87 • class C: WISE sources with 0.10 < Sb < 0.24 and 0.14 < Sg < 0.38. 

88 All the WISE sources with Sb <0.10 or Sg <0.14 are considered outliers of the WGS. Sources of 

89 class A are the rarest with respect to the other classes (Paper IV). 

90 Our association procedure consists in ranking all the WISE sources within the searching region 

91 of an unidentified 7-ray source as described above and indicating as a 7-ray blazar candidate the 

92 positionally closest source belonging to the highest class. Our association procedure provides 

93 a completeness of 87% based on the a posteriori re-association of the ROMA-BZCAT blazars, 

94 detected by WISE and associated with Fermi-LAI sources. 



95 4. 7-ray blazar candidates among the unidentified INTEGRAL sources 

96 We applied our new association method to the case of the 86 UISs selected above. This 

97 process allows us to verify if there is a 7-ray blazar candidate within the positional error region of 

98 each UIS analyzed. 

99 We considered a searching region with radius equal to the position error at 90% confidence 

100 level, as reported in the 4IC catalog; then, we estimated the IR WISE colors for all the sources that 

101 lie within the searching region. 

102 Running our association procedure, we found that 68 out of 86 UISs have only outliers of 

103 the WGS lying in their searching regions, while within the remaining 18 UISs we found 4 7-ray 



104 blazar candidates of class A, 12 of class B and 2 of class C. In Table 1, we present the list of 

105 7-ray blazar candidates found for the 18 UISs together with their IR colors, as well as the s^ and 

106 Sq parameters. We also estimated the probability to find a generic WISE source with the same s 

107 values in 36 random circular regions of the WISE sky having the same radius R of the searching 

108 regions. We found that this is smaller than 10"^. We note that the positional accuracy of the UISs 
is a least order of magnitude better than that of the unidentified 7-ray sources in 2FGL. 
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Summarizing our results, we found 18 WISE 7-ray blazar candidates that could be candi- 
date counterparts of the corresponding UISs responsible for the hard X-ray emission detected by 
INTEGRAL. 



113 5. 7-ray blazar candidates observed by Swift 

114 We found that among the 7-ray blazars selected according to our association procedure, 10 

115 candidates out of 18 have at least one Swift pointed observation. In addition, none of these 7- 

116 ray blazar candidates has a 7-ray counterpart in the 2FGL. We reduced and analyzed these Swift 

117 observations to verify if these WISE 7-ray blazar candidates have an optical-UV or soft X-ray 

118 counterpart. Here we report the data reduction and analysis procedures used in our Swift data 

119 analysis. The comparison between the Swift and the WISE images will be presented in Section |6l 



5.1. UVOT data analysis 

We followed the same UV-Qptical Te lescope (UVOT) reduction procedure described in 



22 (|Tramacere et al. 20071: iMassaro et al. 2008al) consequently, here we report only the basic details. 



123 Several filter combinations are available for UVOT observations; however, we note that not all 

124 the optical and UV data are available for each source. The detection algorithm UVOT DETECT 

125 was used to confirm the presence of the optical-UV counterpart of the 7-ray blazar candidates. We 

126 then performed the photometric analysis using the UV OT SOU RC E tool. Counts were extracted 

127 from a 6" radius aperture in the V, B, and U filters and from a 12" radius aperture for the other 

128 UV filters {UVWl, UVM2, and UVW2), to properly take into account the wider Point Spread 

129 Function in these bandpasses. The count rate was corrected for coincidence loss, and the back- 
ground subtraction was performed by estimating its level in an offset circular region at 20" from 
the source. 



132 The correction for the interstellar reddening was obtained assuming the E(B — V) values 

133 from Schlegel et al. (1998) and the corrections described in Cardelli et al. (1989), while the fluxes 



134 were derived with the same conversion factors given by Giommi et al. (2006). 



5.2. XRT data analysis 

The X-Ray Telescope (XRT) data reduction used in the following is also the same one de 



137 scribed in (|Tramacere et al. 20071 IMassaro et al. 2008bl: iMassaro et al. 2011c|) : here we only re 



138 port the basic details. 

139 The XRT data analysis has been performed with the XRT DAS software (v. 2.1), developed 

140 at the ASI Science Data Center (ASDC) and distributed within the HEAsoft package (v. 6. 10.0). 

141 Event files were calibrated and cleaned with standard filtering criteria using the xrtpipeline 

142 task, combined with the latest calibration files available in the Swift CALDB distributed by 

143 HEASARC. Only events in the energy range 0.3-10 keV with grades 0-12 were used. When 

144 more than a single Swift pointing of each source has been performed and is available within the 

145 Swift archive, we combined several low S/N observations, because the the co-added X-ray im- 

146 age increases significantly the source detection. No signatures of pile-up were found in our XRT 

147 observations. 

148 Given the low exposure of the Swift observations it was not possible to carry out a detailed 

149 spectral analysis, so unless stated otherwise we used the detection algorithm detect, a tool of the 

150 XIMAGE package for all the Swift observations. The detect algorithm locates the X-ray point 

151 sources using a sliding-cell method taking into account the average background intensity. The 

152 position and intensity of each detected source is calculated using a box whose size maximizes the 

153 signal-to-noise ratio. This detection algorithm has bee n extensively used in the Swift serendipitous 



survey in deep XRT gamma-ray burst fields (see also IPuccetti et al. 201 ll for additional details). 



155 Statistical and systematic uncertainties on count rates are added quadratically. 

156 Finally, we measured the net count rates for each detected soft X-ray source and we converted 

157 them into fluxes assuming a power-law spectrum with spectral index 1 and using WEBPIMM^ 



158 6. Results on the Swift analysis 

159 Applying our new association procedure developed for the unidentified 7-ray sources of 

160 Fermi-hAT to the UISs, we found that 18 sources out of the 86 analyzed have a 7-ray blazar 

161 candidate as possible counterparts. We note that this new association method proposed for the 



*http://heasarc. nasa.gov/Tools/w3pimms.html 



162 unidentified 7-ray sources of Ferm/-LAT does not have the same efficiency when applied to soft 

163 7-rays and/or hard X-rays. In fact, in th e FermZ-LAT energ y range (i.e., 30 MeV - 10 GeV) blazars 



are the largest known 7-ray population (|Nolan et al. 2012|) . while the hard X-ray band is generally 



165 dominated by the emission of different classes of AGN, such as Seyfert galaxies, which constitute 

166 ~ 24% in comparison with the 2.4% for blazars already associated in the 4IC. This implies that 

167 the WISE 7-ray blazar candidates are not necessarily the low-energy counterparts of the UISs. 

168 The relation between the IR spectral shape in the WISE energy range and that in the 7-rays 

169 is based on our association method (e.g.. Paper II). On the other hand, in hard X-rays there is not 

170 yet evidence of a correlation between the IR and the X-ray emission of blazars, thus the eventual 

171 association between the WISE 7-ray blazar candidate and the UISs is less robust that in the case of 

172 the Ferm/-LAT sources. 

173 For 10 out of the 18 WISE 7-ray blazar candidates, we also found optical-UV and X-ray ob- 

174 servations available in the Swift archive that could be helpful to verify if they are the low-energy 

175 counterparts of the UISs. We found that 7 out of the 10 WISE 7-ray blazar candidates in the 

176 Swift archive have a clear counterpart in X-rays and in the optical-UV bands, showing a typical 

177 Spectral Energy Distribution (SED) dominated by non-thermal emission, as for the two cases 

178 shown in Figure [T] where J03565 1.52-1-624553. 8 has also a radio counterpart at 1.12 arcsec from 

179 the WISE position. In Table 1, we report the INTEGRAL name together with the WISE 7-ray 

180 blazar candidates, the J2000 coordinates RA and DEC, the distance between the WISE source and 

181 the INTEGRAL position in arcsec, the WISE colors (i.e., c^ = [3.4] - [4.6], C23 = [4.6] - [12], 

182 C23 = [12] — [22]), the S5 and the Sg derived from our WGS method, the Swift UVOT detections and 

183 the Swift XRT detections with the X-ray counts in the soft (0.3-1 keV) and in the hard (i.e., 1-10 

184 keV) bands together with the X-ray hardness ratio HR derived from the net number of counts. 

185 For remaining three WISE 7-ray blazar candidates we did not find a clear counterpart in Swift 

186 observations. This result could be due to the short exposures of the archival observations. 

187 However, we note that in the above 10 candidates, the Swift XRT flux is not in agreement with 

188 the extrapolation of the INTEGRAL spectrum, which is generally one order of magnitude larger 

189 than the Swift XRT estimate. This discrepancy is not sufficient to exclude the blazar association 

190 of the UISs because blazars exhibit rapid X-ray variability; however, the it could suggest that the 

191 blazar is not the most probable low-energy counterpart for the UISs, in agreement with the fact that 

192 they are not the dominant class of AGN in the hard X-rays. We note that the 7-ray blazar candidates 

193 found with our method are WISE sources, detected in all four WISE bands, in particular at 22 /im 

194 as the case of IGRJ14549— 6459 shown in Figure |2l for which the WISE candidate counterpart 

195 appear to have the IR colors of blazars. 

196 In Figure |3] and Figure IH we also show the comparison between the soft X-ray and the optical- 
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Fig. 1. — The SEDs of two example of WISE 7-ray blazar candidates: J03565 1.52+624553. 8 
and J06523 1.40+533 13 1.1. The dashed line in the left panel is the typical log-parabolic 
model adopted to describe the non-thermal SED of J03565 1.52+624553. 8 while in the case of 
J06523 1.40+533 13 1.1 a simple power-law, over 6 orders of magnitude, has been used. As de- 
scribed in Section |6] there is discrepancy between the XRT fluxes and those of INTEGRAL that do 
not support the blazar association of the UIS. 




Fig. 2. — The WISE IR images at 3.4/im, 4.6/xm, 12//m, 22//m, respectively for the Field of View 
(FOV) for IGRJ14549— 6459, centered on the position of the 7-ray blazar candidate. The WISE 
7-ray blazar candidate is highlighted with the green cross in the center of the WISE images. It is 
clear that the source selected with our association method is the only one detected in all 4 WISE 
bands. 



197 UV images of Swift with the WISE IR data, for the FOV of two examples of UISs centered on the 

198 positions of our 7-ray blazar candidates: IGRJ06523+5334 and IGRJ13045— 5630, respectively. 



199 7. Summary 

200 We applied our new association method successfully used for the unidentified 7-ray sources 

201 in the 2FGL to the UISs to test if it is possible to find WISE blazar counterparts at low energies 

202 responsible for the hard X-ray emission detected by INTEGRALwithin the searching regions of the 

203 UISs. 

204 We found that 18 out of 86 UISs analyzed clearly have a blazar counterpart within the search- 

205 ing regions, and for 10 of them also Swift archival observations are available. However, for the 

206 latter 10 sources in Swift we did not find a good agreement between the Swift X-ray flux and the 

207 one estimated by the extrapolation of the INTEGRAL spectrum. Thus, we are not able to confirm 

208 if the WISE 7-ray blazar candidates found with our method could be associated with the UISs 

209 considered. 

210 On the other hand, one crucial result arises from our analysis. We found that the WISE 7-ray 

211 blazar candidates selected from our method in these serendipitous Swift observations of the UIS 

212 fields of view have clear optical and/or UV and soft X-ray counterparts. This is strongly in agree- 

213 ment with the expectations driven by their blazar nature. It is worth noting that J03565 1 .52-1-624553.8 

214 has also a radio counterpart and in addition, the SEDs of these WISE 7-ray blazar candidates are 

215 in agreement with a non-thermal shape over several orders of magnitude. 

216 Finally, we remark that future follow up observations, in particular spectroscopic optical data, 

217 are necessary to clari fy the nature of the WISE 7-ray blazar candidates and conseq uently the nature 



218 



of the UISs (see e.g. JMasetti et al. 2008l : lMasetti et al. 20091 : Masetti et al. 2010h 
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Fig. 3. — The WISE 3.4 /xm IR contours (green) overlaid on the Swift optical-UV and X-ray 
images, for the FOV of IGRJ06523+5334, centered on the position of the selected WISE blazar 
candidate. It is clear that the WISE blazar candidate (red arrow) has a clear counterpart in the soft 
X-rays and in the optical-UV bands. 




Fig. 4.— Same of Figure |3] for the IGRJ 13045-5630 FOV. 



